Logarithmic divergence of the block entanglement entropy for 
the ferromagnetic Heisenberg model 
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Recent studies have shown that logarithmic divergence of entanglement entropy as function of 
size of a subsystem is a signature of criticality in quantum models. We demonstrate that the ground 
state entanglement entropy of n sites for ferromagnetic Heisenberg spin- 1/2 chain of the length L in 



a sector with fixed magnetization y per site grows as i log. 



n(L — n) 



C(y), where C(y) = 2ne^ - y 2 ) 



PACS numbers: 

Recently it has been argued, on the example of the 
exactly solvable antifcrromagnetic Heisenberg spin 1/2 
chain 



H 



xxz 



= J 



E 

(=1 



Aafaf 



(1) 



that for critical (gapless) quantum system (for the XXZ 
model when A belongs to the interval (—1,1)) the en- 
tanglement entropy of a block of n spins diverges loga- 
rithmically as 7 log 2 n, while for non critical systems (A 
outside the above mentioned interval), it converges to 
a constant finite value 0, 0, 13 • This property was in- 
terpreted in the framework of conformal field theory |J| 
associated with the corresponding quantum phase tran- 
sition and the prefactor 7 of the logarithm related to the 
central charge of the theory c = 37 (for the XXZ model 
this gives 7 = 1/3). 

The aim of this Letter is to show that the entanglement 
entropy of a block of spins in the ground state of the 
antiferromagnetic XXZ model at the point A = — 1 
grows faster than for other critical points — 1 < A < 1, 
namely as 7 log 2 n with the logarithmic prefactor 5 < 
7 < 1. 

Our approach uses the permutational invariance of the 
ground state of JIJ at A = — 1, this allowing to com- 
pute the entanglement entropy exactly for blocks of ar- 
bitrary size and system of arbitrary length. To this 
regard we remark that by performing the transforma- 
tion which overturns each second spin along the chain 
(we assume the length of the chain even) the Hamilto- 
nian <£Q) for A = — 1 reduces to the isotropic Heisen- 
berg ferromagnet J2J ■ Since this transformation does not 
change the entropy of entanglement, one can compute the 
block entropy of the antifcrromagnetic Heisenberg chain 
at A = — 1 directly from the one of the isotropic fer- 
romagnetic model. It is worth noting that, in contrast 
with critical points — 1 < A < 1, the point A = — 1 can- 
not be studied by means of conformal field theory since 
this point is not conformal invariant Q , the ground state 



being infinitely degenerated at A = — 1. 

The paper is organized as follows. After introducing 
the model we formulate a theorem which gives the analyt- 
ical expression of the eigenvalues of the reduced density 
matrix. Using this theorem we compute the entangle- 
ment entropy of a block of size n in the finite system of 
total length L for two specific choices of the ground state 
sector. Taking the limit of large subsystem sizes, we de- 
rive analytical expressions for the entanglement entropy 
Sf n ) of a block of spins of size n in the ferromagnetic 
ground state, both for n, L ^> 1 and for n 3> 1, L = 00. 
As a result, we obtain that in the ground state sector 
with a fixed value of S z the block entanglement entropy 
grows for large n, as SV n ) = \ log 2 n \ L ~ n ) ; while in the 
ground state sector in which all the S z components of the 
spin multiplet are equally weighted, S7 n ) = log 2 (n + 1) 
for arbitrary n and L. 

We consider the ferromagnetic Heisenberg model with 
nearest neighbor interaction, 



H 



XXX 



= -J 



E ( cr « cr *+i 



3/ 



(2) 



where a are Pauli matrices, J > denotes the exchange 
constant and L the number of spins (we assume periodic 
boundary conditions L + 1 = 1). As is well known, the 
ground state of belongs to a multiplet of total spin 



S = -§ and is degenerate with respect to S 
1 ^ 
down, 



L _L 
2 ' 2 " r 

In the sector with a fixed number N of spins 

L 

2 ■ 



i.e. with a fixed S z = N — 4, the ground state is 



+ 



(3) 



obtained by the action of the rising operator S + 
on the vacuum state with all spins down 

1*^(^111.4). 

All eigenfunctions J2J correspond to the same ground 
state energy E = of the XXX model 0. The struc- 
ture of the state © is given by 

l^(L,iv)) = ^L=^irT...T4i...i) (4) 

p 



N 



L-N 



2 



where the sum is taken over all possible distributions of 
N spins on L sites and the binomial coefficient = 
n\(l~n)\ takes care of the normalization. Note that |@J 
is also a ground state for the model of interacting bosons 
0, while for the partially asymmetric exclusion process 
ASEP with N particles hopping with hard-core exclu- 
sion on a closed chain of the length L, (@J represents a 
steady-state vector. We will be interested in the ground 
state entanglement (von Neumann) entropy 5(„) of a 
block of n (not necessarily contiguous) spins 



S( n ) = -tr(p n log 2 p n ) = - y] A fc log 2 Afe, 



(5) 



where p n is the reduced density matrix of the block, ob- 
tained from the density matrix p of the whole system by 
tracing out external degrees of freedom p^ — tr^ L _ n jp 
(notice that due to the permutational symmetry of the 
ground state Sr n \ does not depend on the particular 
choice of the block but only on its size n) . In Eq. (JSJ) Afe 
are the eigenvalues of the reduced density matrix which 
are all real, nonnegative, and sum up to one: Y] Afe = 1. 

The density matrix p for a degenerate ground state is 
given by 

L 

p=Y / a N \*(L,N))(*(L,N)\, ^^ = 1, (6) 

N=0 

where ao, oci, is a set of nonnegative coefficients. 

Denoting the reduced density matrix in a fixed sector 
with N spins up by p n (N), 

p n (N) = tr {L _ n} \y { L,N))(y(L,N)\, (7) 

where \*&(L, N)) is given by Q), one can write the gen- 
eral reduced density matrix as 



Pn =^ a N Pn(N). 
N=0 



(8) 



In the following we consider two choices for the coeffi- 
cients {cti}: 



(a) a t = 5in, 

(b) ao = ai - 



L+l 



(9) 
(10) 



(the analysis for arbitrary {an} proceeds in similar man- 
ner). The choice (a) corresponds to the case when a 
small anisotropy single out a sector with TV spins up re- 
sulting in a pure state of a global system, see il ll(il) . The 
choice (b) corresponds to an equilibrated density matrix 
(i.e. with all components of the ground state multiplct 
equally weighted) which preserves the SU (2) invariancc 
of the Hamiltonian J5J) (this case is equivalent to infi- 
nite temperature). Using the general property of the 



entropy of composite systems: S, 



(n) 



(L- 



\, and its 



invariance with respect to the inversion of all spins, we 



can restrict the analysis, without loosing generality, to 
the case: n < -j, N < -j. The computation of the block 
entanglement entropy is drastically simplified by the fol- 
lowing 

Theorem: The eigenvalues of the reduced density matrix 
p n {N) of a block of n spins in the sector with N spins 
up in the ground state of the ferromagnetic Heisenberg 
model 0) are given by 



X k (L,n,N) 



C2C 



N-k 



0,l,...min(n,/V). (11) 



The proof of the theorem follows from the decomposition 
of p n (N) with respect to the symmetric orthogonal sub- 
spaces of the system of n spins, classified by the integer 
k = 0, 1, ... min(n, N) giving the number of spins up in 
the block 



Pn(N) 



min{n, N} 

E 

fc=0 



c k \4>{n 7 k))(ip{n,k)\. 



(12) 



Here \ip{n, k)) denotes the symmetric state with k spins 
up among n spins 



(13) 



i — k 



and Cfe is the corresponding probability = " L k (no- 

tice that C^ZJl is the number of states with k spin up 
in the block of n spins and is the total number of 
states). Expression (|T2"|) can be rewritten as 



min{n,iV} 

i (AT) = Yl x kPn{k) 



(14) 



k=0 



where p n {k) is the density matrix of the state \ip(n,k)) 

C n C Ij ~ n 

and the coefficients Afc = k „L k sum up to one, ^2 X k = 

N 

1. From this it follows that p n (N) is the density matrix 
associated with the ensemble of orthogonal pure states 
{Afc, /?„(£;)} and therefore it has a block diagonal form, 
each block having only one nonzero eigenvalue Afc which 
coincides with the expression Qllfl. This concludes the 
proof of the Theorem. 

We remark that the specific case N = n = ^ was also 
considered in Ref. 8J. Having found the eigenvalues of 
p n (N) one can easily compute the entanglement entropy 
S7 n ) for arbitrary L,n and N. 

Case (a) To obtain an analytical expression for S( n ) , 
from the exact expression (JSJ lllfl . we observe that for 
blocks of large size, n»l, the dominant contribution to 
the sum (jSJ) comes from the eigenvalues Afc with large k. 
In this case one can approximate the binomial coefficients 
in 1|11|) by the normal distribution, see e.g. [j|: 



Gm m n—rn 
n P Q 



1 



y/2iTnpq 



exp 



(m — npY 
2npq 



npq > 1, 
(15) 
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Figure 1: Entanglement entropy as obtained from the ex- 
act expressions Eqs. 15I11H . as a function of the block size 
n and for L — 20 (up triangles), 50 (down triangles), 100 
(squares), 150 (diamonds), 200 (circles). Filled (empty) sym- 
bols correspond to p — 1/10 (p = 1/2). Continuous curves 
represent the analytical expression in Eq. (lib! . For n > L/2, 
S(n) = 5 , ( L _ n) (not shown). 
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Figure 2: Entanglement entropy as function of a number of 
sites involved, for different values of p = yjj, |. Compari- 
son of exact formula (points) with the limiting expression I18H 
(continuous curves). 



where < p < l,q = 1 — p . Using this approxima- 
tion, and defining p = N/L, the eigenvalues can be 
written as 



X k (L,n,N) 



/~ln m k „n — k /~iL — n N —k „L~n—N '+k 
U kP 1 ^N-kP 1 



(JL „N nL — N 



1 1 



: exp 



nP T 



2a 



where a = pq ^ L ■ Substituting this expression into (0 
and replacing the sum with an integral, we obtain 

J R.(log 2 ^j dx, 
1 ( (x-p) 2 ' 



R 



■. exp 



2a 



For large n the limits of the integral can be extended to 
include the whole real axis, after which the result of the 
integration gives 

S [n ){p) « \ log 2 (27re M ) + ~ log 2 < L ~ n \ ( 16 ) 

Notice that this approximate result is valid for npq S> 1 
and in the limit npq — > oo it becomes exact. From 
the analytical expression i|16|) the following properties 
can be easily derived: i) S( n )(p) = S7„)(l — p), ii) 



S(n)(p) = S(i_„)(p), Hi) dS {n) {p)/dn = only at 
n = f , vi) dS( n) {p)/dp = only at p = §, S(n)(p) 
is a monotonically increasing function of the total length 
L. In Fig. Q we compare the exact entropy of finite sys- 
tems, as computed from exact expressions Eqs. I p?) 1111 . 
with the analytical expression (|16Jl . from which we see 
that there is an excellent agreement also for small values 
of npq. In the thermodynamic limit L — + oo, j- — ► p the 



eigenvalues (jl 1|> reduce to 

A fc = CV,CV 



<?>• 



(17) 



and the corresponding entanglement entropy is obtained 
from (|16[) as 

s (n) (p) ~ g lo 8'2 i^epq) + - log 2 n. (18) 

In Fig. 0we plot the exact entanglement entropy of a 
block of size 1 < n < 1000 in an infinite chain JSJ), |T7J|, 
versus the limiting expression l|18fl for different filling p. 
We see that the analytic formula (|18fl gives a good ap- 
proximation even for small finite number of sites n in the 
block. For very small p the convergence is slower (see the 
lowest graph in Fig. |2J) because the validity of formula 
(|15fl crucially depends on the value of npq. 

Thus, for case (o) we conclude that the block entan- 
glement entropy of the ferromagnetic ground state grows 
logarithmically with n, as for critical quantum systems, 



4 



but with a different prefactor, i.e., as \ log 2 n rather than 
\ log 2 n predicted in Q . 

Case (b) In this case the eigenvalues of the reduced 
density matrix are given by 

~k L — k /~tN — n+k 

A fe =-^- y ° L -" = — , fc = 0,l,....n 

N=n-k L 

(19) 

and are independent on k and on the size of the system 
L. The entanglement entropy is obtained as 

S (n) =log 2 (n + l), n = l,2,....L. (20) 

Equations (|11I16I18I ) and i|19l2U|) , corresponding to the 
cases (a) and (b) considered above are the main results 
of the paper. 

It is worth to note that, due to the permutational in- 
variance of the ground state, for any choice of the density 
matrix JHJ the reduced density matrix for a block of size 
n has exactly n+1 nonzero eigenvalues (see the theorem) 
in the ground state. This implies the upper bound for the 
entropy S max (n) = log 2 (n + 1), which is achieved in the 
case of a thermally equilibrated density matrix (case (b)). 
The lower bound of logarithmic growth Sr n ) ~ \ l°g2 n 1S 
achieved for the "anysotropic" choice corresponding to a 
pure state © of the whole system (case (a)). For generic 
choice of the coefficients {aAr} in Kil.SI . S/ n \ will grow as 
7 log 2 n with ^ < 7 < 1 . 

We also note that i|20|) is a monotonically increasing 
function of n, attaining maximum for the whole system 
n = L, while in the case of pure state the maximum is 
achieved for a block of half-system size n = \. This 
feature is related to the fact that the ground state of a 
ferromagnet is highly degenerate and the total system for 
the choice (|l()fl is in the maximally mixed state. 

Another remark concerns the origin of the logarithmic 
prefactor 7 = ^ in formula l|18|l . Apparently 7 is not 
related to any central charge in since A = — 1 is not a 



conformal point. We find that in our case the prefactor 
7 is related to the spin s per site, i.e., one can show 
that for a ferromagnetic spin s chain (i.e. with on-site 
spin s), the block entanglement entropy in the ground 
state sector grows like ~ const + s log 2 n (details 
will be presented elsewhere) . We finally remark that it is 
of interest to generalize Eas. (|16l20l) to the case of nonzero 
temperature, where excited states have to be taken into 
account. Work in this direction is in progress. 
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